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ABSTRACT 



A compact heat exchanger testing facility was constructed at the 
United States Naval Postgraduate School. The testing facility is used 
to collect heat transfer and friction factor data by the steady state 
steam to air technique. Heat exchangers with frontal dimensions of up 
to 12 inches by 12 inches can be tested. The maximum pressure drop which 
can be produced across the air side of a heat exchanger is 2.5 psi. The 
testing facility was designed to obtain data for exchanger air side 
Reynolds numbers from 500 to 10,000. A description of the apparatus, the 
experimental method and equations and an evaluation are given. 
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- Test core flow length, ft. 

- \/2 h/k(T 

- Mass rate of air flow, Ibjjj/hr. 

- Pressure, in. H 2 O, in. Hg, Ibf/in^, Ib^/ft^ 

- Barometric pressure, in. Hg 

- Heat transfer rate, BTU/hr 

- Gas constant for air, 53.3 ft lbf/(lbjjj °R) 

- Hydraulic radius, (A^|,/A^) , ft, (4 rj^-hydraulic diameter) 

- Mass rate of steam flow, lb „/hr 
Temperature, °F 

- Absolute temperature, °R, ok 

- Unit overall thermal conductance, BTU/(hr °F ft^) 

- Velocity, ft/sec 
Specific volume, ft^/lb^ 

- Mass rate of condensate from test core, lb|jj^/hr 
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X 



Ratio of pressure upstream of the fluid metering orifice to 
pressure differential across the orifice 



Xc 
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Humidity correction to the specific heat of 


air 


Xm 
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Humidity correction to the density of air 
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Net expansion factor for a square-edged metering orifice 
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Letter Symbols 
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- 


Ratio of fluid metering orifice diameter to 


duct diameter 
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Denotes difference 
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- 


Fin thickness, in. 
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Fin temperature effectiveness 




^0 


- 


Total surface temperature effectiveness for 


air side 


^s 
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Total surface temperature effectiveness for 


steam side 


cr 
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Ratio of free-flow to frontal area, 
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Dynamic viscosity, Ibm/ (hr ft) 
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Density, lb 
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Average density of air in test core for isothermal test, 
lb ^/ft3 


^m 
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Mean density of air in test core for hot core test, Ib-'^j^/ft 


Dimensionless Groupings 




NRe 
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Reynolds number, (4r^G^) 
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Stanton number, (h/GCp) 




Npr 
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Prandtl number, (^Cp/k) 
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- 


Colburn- j (Ng^ Np^2/3) 




NTU 
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Subscripts 



1 


- Air upstream of core 


2 


- Air downstream of core 


o 


- Air metering orifice 


a 


- Air side of test core 


c 


- Core 


f 


Fin 


f 


- Air film in test core 


s 


Steam side of test core 


so 


Steam metering orifice 


si 


Steam in top header 


s2 


- Steam in bottom header 


w 


Wall between air and steam sides of test core 


s 


Square fin 


t 


- Triangular fin 


c 


- Circular tube 



Force and Mass Units 

Ib^ - Denotes pounds mass in distinction to 



ibf 



Denoting pounds force 
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Introduction 



Steady state tests transferring heat from condensing steam to 
flowing air are used to determine experimentally the heat transfer 
characteristics for compact heat exchangers. The friction factor is 
determined from the air pressure drop across the heat exchanger. Steam 
and air are used as the working fluids in the tests because they are 
readily available and produce effects which are easily analyzed. The 
heat transfer characteristics are given by Colburn-j heat transfer 
modulus evaluated from the steady state steam to air test (referred to 
here-in as hot core test). The friction factor is evaluated from both 
the hot core test and the isothermal test (referred to as cold core test). 
The Colburn-j heat transfer modulus and friction factor are normally 
presented graphically as functions of the heat exchanger Reynolds number. 
Development of these three parameters is given in the theory section. 

In addition to determining the Colburn-j and friction factor, a 
heat balance is performed on the exchanger to check the consistency of 
the experimental data. 

The testing facility is similar in design to one built at Stanford 
University in 1947. The Stanford facility was built under a United 
States Navy contract to evaluate finned surfaces for use in compact hdat 
exchangers and is described by Kays in reference [3]. 

Compact heat exchangers which are suitable for steady state steam 
to air testing have plate-fin geometry as shown below- 
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The steady state steam to air testing facility was designed and 
built to accommodate a heat exchanger with frontal dimensions of up to 
12 inches by 12 inches. The facility is capable of producing air side 
Reynolds numbers from 500 to 10,000. An air side pressure drop of 2.5 
psi can be obtained. 

The testing facility is described and the experimental method and 
equations for data reduction are presented. The data reduction equations 
are conveniently arranged for programing on a digital computer. A test 
of a heat exchanger of known characteristics is presented as an evalua- 
tion of the testing facility. The operating procedures and the computer 
program for data reduction with accompanying instructions are given in 
the appendix. 
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2 . Theory 



Heat Transfer Characteristics . The heat transfer characteristics 
of a compact heat exchanger are best represented by dimensionless group- 
ings. This generalizes the results so that they may be used for working 
fluids other than steam and air. The succeeding development is that 
commonly used in the field; see Kays and London [4]. The Colburn-j heat 
transfer modulus is selected for this purpose as it combines heat transfer 
characteristics and fluid properties. The Colburn-j modulus is defined as 




where the ,unit conductance for convective heat transfer, h, must be deter- 
mined experimentally. The use of a steam to air system is one of the 
most direct methods of experimentally determining h. 

The unit conductance for convective heat transfer, h, for the air 
side, air being the fluid of interest, is determined from energy balance 
considerations within the exchanger. Low pressure slightly superheated 
steam, about 6 psig with 5 degrees superheat, enters the top of the heat 
exchanger steam side. The steam reaches its saturation temperature in a 
very short distance after entering the exchanger. Constant temperature 
is maintained as the steam condenses for the remainder of travel. The 
air is heated continuously as it passes through the exchanger. The 
sketch below illustrates the temperature variations in the working fluids 
as a function of the heat transfer area. 




An energy balance for the differential area dA yields the equation 



for the heat transfer rate as: 



dq 

(heat transfer 
rate) 



ft Cp dXair = 

(heat transfer rate 
to the air) 



U d A AT 

(heat transfer rate 
from steam to air) 



Rearrangement gives: 



dTair 



(Tsteam-Tair) 



= _ U d A 

m c 



For small over-all changes in fluid properties and a relatively uniform flow 
distribution in the heat exchanger, the over-all unit conductance, U, 
may be assumed to be constant. Integration over the total area then 
yields : 



in 



(- 



Tsteam - Tair 



ir out ^ _ U 

ir in y m 



Tsteam - Tair 

The dimensionless grouping on the right side of the equation is defined 
as the Number of heat Transfer Units, NTU. Thus the equation may be 
written as: 



in 



(- 



Tsteam - Tair 



Tsteam - Tair 



ir out \ 
ir in J 



= - NTU 



or 



(- 



Tsteam - Tair 



ir out \ _ 
ir in / 



-NTU 



Tsteam - Tair 

The over-all unit conductance based upon the total air side heat 
transfer area must include a factor to account for the temperature 
variation in the fin areas. The over-all temperature effectiveness, n, 

may be written as: ^ 

n = 1 - (1 -nf) 

A 
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where Af is the fin area and7[f is the fin temperature effectiveness. 



The effectiveness for a fin of constant conduction cross section as 



given by Kays and London [4] is: 



tanh 






where 




Z = one half fin height for wall to wall fins 



Summing the resistances to heat flow then yields: 



1 

U 



^ A^ ?a*i. 

Ag hg A„ ki. 



1 




(over-all (steam side (wall 

resistance) resistance) resistance) 



(air side 
resistance) 



The convective coefficient for the air side, h, may be solved for giving: 



A value of hg equal to 2000 BTU/(hr ft^ °F) is assumed. A fifty 
percent error in this assumption results in approximately two percent 
error in the value of h. Accordingly the heat transfer coefficient for 
the air side, h, may be evaluated by experimentally measuring the fluid 
temperatures and the air mass rate of flow. 

Friction Factor . The flow of a fluid through an exchanger results 
in a pressure differential caused by the flow contraction as the entrance, 
the fluid friction in the exchanger and the flow expansion at the exchanger 
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h 




( 1 ) 
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exit. These pressure effects are graphically represented below. 



HEAT 

1 EXCHANGER 2 




To define a meaningful heat exchanger friction factor, the entrance 
and exit pressure effects must be elimated from the total pressure 
differential. The friction factor will then be based solely upon the 
mean wall shear. The problem of entrance and exit pressure effects has 
been investigated both analytically and experimentally by Kays. [2]. 

The effect of pressure loss in the downstream duct is included in 
the total pressure drop in order that point 2 may be located downstream 
sufficiently to assure complete pressure recovery. The total pressure 
differential from point 1 to point 2 is then: 

^total “ ^^a ^^c " ^^b ^^cluct 

The expression for the pressure effect due to flow contraction at 
the entrance is given by Kays [2] is: 

A Pa ^ + Kcv| 

2gc', 2 9c 
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and for the flow expansion at the exit is: 

^__^(l-<r 2 ) _ KeVg 

?b 

with denoting the exchanger free flow to frontal area ratio. ind 

Kg are contraction and expansion coefficients which are functions of 
exchanger geometry and Reynolds number. Noting that the velocity may be 
related to the mass velocity and that the density change from section one 
to a and b to two will be very small compared to the total density change, 
the above equations yield: 



and 



Ap.,= „ (|-(t 2 - Ke^ 



The pressure drop in the exchanger may be found by writing Bernoulli's 
equation in the differential form to include the Fanning friction term. 
Using the mass velocity, G, the equation as given by Kays [3] is: 

dPr \ \ . - / dx 



or 



— C — 

? 29 c 



-dPc = -S 






+ f 



2gce^ 



2g< 






, 



By defining the mean specific volume in the exchanger as: 

•L 

- I±.\ - ± ' 






.ml 



I 



and integrating from p.^, 0 to L, the exchanger pressure differential 
is given by: 



AP - 

^ 2 gcea 



Pc 



Pb 
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The mean specific volume is found from its definition by replacing it 



with its value as given in the equation of state for a perfect gas: 

P 



j_ ^ 



RT 



and expressing temperature and pressure as functions of x. As the 
pressure variation is small compared to the temperature variatioi^:, 
arithmetic average value of P may be used: 

p = P-a + 

1 

The expression for temperature as a function of x may be written from 
the temperature-NTU relation as: 



NTU 



X 



T = Tg - (T^ - T^) e 



where T^ is the steam temperature. Substitution of the expressions for 
temperature and pressure into the equation for the mean specific volume 
yields: 






_L 2R 

em= L(P^ + P ) 



0> 



NTU 



Ts-(Ts-Ta)e‘ 



dx 



Integration and using the equation of state for a perfect gas to simplify, 
and assuming no temperature change from one to a and b to two the follow- 
ing equation obtains: 



-L = 2 

Pm +"TbPb 






(izzjd 



NTU 



An alternate expression may be written by substituting the NTU-tempera- 
ture relation into the above equation. This expression is: 

J . 2R (Tj -LMTd) 

?m” Pd+Pb 
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with LMTD being the Log Mean Temperature Difference. 



The pressure loss in the downstream duct may be evaluated as: 

^Pgiuc t- - 4 \ Md 

\ Dd / 2gc 

where the density, p^, the friction factor, f^, the length of duct 
from section b to 2, the hydraulic diameter, D^, and the velocity, 

are all quantities associated with the downstream duct. The velo- 
city in the downstream duct may be related to quantities appearing in 
the previous pressure loss equations. If the downstream duct has the 
same dimensions as the exchanger and small density variations are 
neglected, the continuity equation gives 



m 



lAcjlfbllAfn 



or 



Vd = 



G 0~ 
?b 



Thence 






@2 0-2 



Dd / 2gceb 



Again assuming that the density changes between section 1 to a and 



ap = ® 



2gc 



b to 2 are negligible, the various pressure drops are summed yielding: 

'' (|-(r 2 + Kc) , 

- e, ^ e, e2 

e; +AEZ_ 



Manipulation of the equation gives: 






2gc 



f L 



he, 



-h 



m 



Kc -(o-2-l- l) 
— 



-l- 



(j + Kg-f- 0-2 ( I + 4 D<^' ) 



and so 



Iving the equation for the friction factor we obtain: 



f = 



_ r~h 
L 



^toi-al _ (Kq- CT^ -I) 



[l+Ke+ <r2(l+4p^ 1" 




?2 J 





( 2 ) 



Therefore the friction factor, f, for a heat exchanger of given geometry 
may be found by experimentally determining the mass rate of air flow, 
the entrance and exit air pressures and temperatures and the steam 
temperature. 

Reynolds Number . For uniformity of reporting test results, the 
Reynolds number is defined in the same manner as given by Kays and 
London. [4]. This definition embodies the hydraulic diameter as the 
characteristic length and the velocity based on the minimum free flow 
area. Thus: 



No. = ^ 






Ac 



To evaluate the foregoing characteristics, the test apparatus 
must be capable of accurately measuring the air state upstream and 
downstream of the exchanger, the steam state within the exchanger and 
the mass rate of air flow. 



3. Description of the Test Apparatus 

General . Air is induced through one side of the compact heat 
exchanger by a centrifugal compressor and associated ducting. Steam 
is passed through the other side of the exchanger by the steam system. 

The air and steam systems are shown diagramatically in Fig. 4 and Fig. 5. 
The steam system is designed to give close control of the steam state 
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within the core. Pressure and temperature instrumentation are pro- 
vided for measuring the air properties upstream and downstream of the 
core, the mass rate of air flow, the condition of the steam entering 
and leaving the core and the mass rate of excess steam and condensate 
leaving the core. 

The test apparatus is shown pictorially in Figs. 1, 2, and 3. 

Air System . The air system is made up of 16 gage galvanized 
steel ducts joined with 1/4 inch steel flanges. The ducting prior to 
the test section has a maximum dimension of 12 inches square and the 
ducting after the test section is 14 inches in diameter. The lengths 
of ducting were designed in accordance with the ASME Power Test Code 
[7] for establishment of uniform flow and are shown in Fig. 4. 

The entrance to the test duct is on the exterior of the building 
to minimize thermal gradients in the air stream. The entrance is three 
feet square and is covered by a fine mesh screen to eliminate foreign 
matter from the air stream. The transition from three feet square to 
12 inches square is made in the entrance section, each side of which has 
the curvature of a quarter ellipse. After a five foot length of duct, 
a second transition is made from 12 inches square to the test core 
dimensions. This is followed by the forward instrument section which 
contains four thermocouple taps, a piezometer ring and two pitot tube 
taps for conducting vertical and horizontal velocity surveys. Follow- 
ing the test core is the after instrument section which contains nine 
thermocouple taps and a piezometer ring. A transition is then made to 
a 14 inch outside diameter duct at the end of which is a standard ASME 
square edged orifice. Piezometer rings are located one diameter 
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upstream of the orifice and one-half diameter downstream of the orifice. 
The duct section downstream of the orifice also contains a thermocouple 
tap. 

Air flow is regulated by a manually operated double sliding plate 
valve located at the turbo-compressor inlet. The valve plates slide on 
Teflon runners and the drive screw is ball bearing mounted to eliminate 
maintenance problems. The double plates constrict the flow symmetrically 
and facilitate fine control of the air flow. Coarse control of the air 
flow is accomplished with the blast gate (discharge valve) of the cen- 
trifugal compressor. 

Air is induced through the test section ducting by a centrifugal 
compressor rated 6000 cfm at 2.5 psi back pressure. The compressor is 
driven by a three phase, 220 volt, 100 hp motor. The compressor dis- 
charges air through a 20 inch square duct to the exterior of the 
building. 

The air flow is metered by means of a standard ASME square edged 
orifice designed to the specifications given in the ASME Power Test 
Code [7]. Four orifice plates made of 1/4 inch type 304 stainless steel 
were designed with the following orifice diameters. The preceeding duct 
has an inside diameter, D, of 13.875 inches. 



d (in) 


R - d 
^ D 


10.406 


0.75 


6.244 


0.45 


3.468 


0.25 


2.081 


0.15 



These diameters were selected so that flow rates of from 250 lbs per hour 
to 22,000 lbs per hour can be measured with the ranges overlapping. The 
orifice plates are fitted with a centering grove on the bottom and a 
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scribe mark on the top so that they may be accurately positioned. 

The ducting downstream of the test section is insulated with two 
inch fiber glass insulation covered by aluminum foil, and the air duct 
is insulhted from the test core by 1/8 inch Teflon gaskets. 

Steam System . Heating of the test core is accomplished by passing 
saturated steam through one side of the core. Before entering the core 
the steam is reduced to a low pressure, about six psig, and a slightly 
superheated state, about five degrees superheat. In a superheated 
condition, the state of the steam can accurately be established by 
pressure and temperature. The few degrees of superheat are quickly 
lost in the test core and have negligible effect on the calculations. 

A considerable excess of steam is passed through the core to prevent 
a thick film of cpndensate from building up on the heat transfer sur- 
faces . 

Steam is supplied from the central steam heating system at 125 
psig in a saturated state. It is strained and then reduced to about 
30 psig through a 1 - 1/4 inch air operated, pilot controlled pressure 
reducer. At this point water is injected, if necessary, to desuperheat 
the steam and then it flows through a 1 - 1/4 inch centrifugal separator 
which removes all entrained moisture. The steam then passes through a 
second 1 “1/4 inch air operated, pilot controlled pressure reducer 
which reduces it to about six psig. This second reduction slightly 
superheats the steam. Just prior to entry into the core the steam 
passes through a short transition section. To provide uniformity of 
flow in the test core, the transition section is filled with aluminum 
shavings held in place by a fine mesh stainless steel screen. Steam and 
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condensate from the bottom of the test core flow through a second 
transition section to a 1 - 1/4 inch centrifugal separator where the 
condensate is removed. A float type trap receives the condensate and 
discharges it at a realitively constant flow rate. The cbndensate is 
then subcooled in a small tap water counter flow heat exchanger and 
collected in a bucket for weighing. The exhaust steam flows from the 
separator through a standard ASME square edged orifice with flange 
pressure taps and is discharged to the atmosphere. 

Two orifice plates were selected to measure the steam flow rates 
and have the following orifice diameters. The inside pipe diameter, 
Ds, is equal to 1.25 inches. 

ds 

ds (in) = dJ 

0.700 .560 

0.890 .712 

The metering orifice was installed as specified in the ASME Power 
Test Code [ 7 }. 

Pressure and thermocouple taps were installed in both the upper 
and lower transition sections and a thermocouple tap was installed 
upstream of the steam orifice. The entire steam system is well insu- 
lated to minimize heat loss. 

Compressed air for actuation of the two steam pressure reducers 
is supplied by two air pressure regulators located on the instrument 
panel. This combination is designed to hold the downstream steam 
pressure to 1 0.2 psig. An air line ftom the air supply of the first 
pressure reducer was fitted into the steam strainer clean-out plug so 
that the entire steam system can be blown dry upon completion of a 
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series of test runs. This should minimize rusting when the system stands 
idle for long periods of time. 
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Pressure Instrumentation . Pressure instrumentation is provided 



for measuring , the gage pressure in the air duct upstream of the test 
core and upstream of the air orifice, the air pressure differential 
across the test core and air orifice, the steam header pressures and 
the pressure differential across the steam orifice. Air pressures are 
measured with well-type single leg water manometers and inclined draft 
gages. Steam pressures are measured with well type single leg mercury 
manometers. 

A one inch and three inch inclined draft gage and a 60 inch water 
manometer are used for measuring the air pressure differential across 
the test core. A 30 inch water manometer is used for measuring gage 
pressure upstream of the test core. The pressure tap downstream of the 
test core is located so that full pressure recovery of the air stream 
is attained. The air orifice pressure differential is measured with a 
three inch inclined draft gage and a 30 inch water manometer. A 60 inch 
water manometer is used to measure the gage pressure upstream of the 
air orifice. Each pressure tap on the air duct consists of four 1/16 
inch holes drilled symmetrically around the duct connected together by 
a piezometer ring of 1/4 inch copper tubing. One quarter inch plastic 
tubing is used to connect the piezometer rings to the manometers and 
draft gages and each manometer and draft gage is isolated by a brass 
toggle valve. 

Three 30 inch mercury manometers are used to measure the gage 
pressure of the top and bottom steam headers and the pressure differential 
across the steam orifice. Water is used as the intermediate fluid 
between the steam and mercury and a constant head of water is maintained 
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by using water potSo The change of level in the water pots is the same 
order of magnitude as the change in level in the mercury wells and these 
changes have insignificant effect on the pressure readings* The water 
pots are stainless steel cylinders three inches in diameter by three 
inches high and are located above the top of the manometers so that air 
pockets will not form in the water lines connecting them to the manometers. 
The connecting water lines are of 1/4 inch stainless steel tubing and 
each manometer is isolated by a stainless steel toggle valve<, The water 
pots are connected to the steam headers and steam orifice by 5/8 inch 
copper tubing. The large diameter tubing was selected so that conden- 
sate can flow freely back into the steam system and not form pockets 
which would change the static head on the manometers. 

Temperature Instrumentation . All temperatures are measured with 
copper-constantan thermocouples. The thermocouples are in assemblies 
manufactured by Honeywell under the trade name MegapaK. The MegapaK 
consists of a **head” for connepting the extension wires, a "sheath” of 
1/8 inch stainless steel tubing with the insulated leads extending down 
through the center and the measuring junction at the end of the sheath. 

The thermocouples used for measuring the air stream temperatures are 
"exposed," meaning that the junction is extended one sheath diameter 
past the end of the sheath. Those used for measuring the steam tempera- 
tures are "remote," meaning that the j Unction is one sheath diameter 
from the end of the sheath and the end of the sheat is sealed for pressure 
tightness. All extension wire is polyvinyl covered 24 gage copper-con- 
stantan. 

The thermocouple voltages are measured with a Honeyx^ell 24 point 
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Multipoint recorder. The reference junction is an ice junction and is 
included in the measuring circuit. The recorder is a D.C. millivolt 
recorder with a scale of zero to ten millivolts. The smallest division 
on the chart paper is 1/10 of a millivolt and temperatures may be 
estimated to 1/100 of a millivolt, The recorder has several modes of 
operation. On the Multipoint setting it will balance and print each of 
the 24 points consecutively at a rate of about two seconds per point. 

On the Selecto-Print setting the recorder will balance and print at 
about two second intervals any number of the 24 points. On the Hold-on 
setting it continuously balances on any one of the 24 points but will 
not print. This setting may be used for obtaining a pen trace by 
attaching a clip-on pen to the print mechanism. The recorder should 
be adjusted for zero and span before use and the adjustment procedure 
is given in the operating procedures. 

The temperatures of the air upstream and downstream of the core, 
and the downstream of the orifice are recorded. Steam temperature in 
the top and bottom headers and upstream of the steam orifice are also 
recorded. To obtain an average temperature of air upstream and down- 
stream of the core a system of thermocouples is used. Upstream of the 
core this consists of three thermocouples placed on the horizontal center 
line. Downstream of the core nine thermocouples are arranged in a three 
by three grid with equal spacing between each of the thermocouples and 
the thermocouples and the wall. In this manner small variations in 
temperature across the stream can be averaged out. All other tempera- 
tures are measured by a single thermocouple. Each penetration through 
the duct for a thermocouple is fitted with a Swagelok compression fitting 
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so that the penetration length can be adjusted. 

The thermocouples are numbered as shown in Figs. 4 and 5 and the 
same number is printed beside the temperature reading on the multipoint 
recorder chart. 



4. Experimental Method and Equations 

General . The experimental method and equations used are essentially 
those outlined by Kays [3] for a similar facility. These equations are 
suited for use with tables and graphs. The equations outlined herein are 
suited for use in a digital computer program which will accept raw test 
data and produce finished results. The equations used in preparation of 
original graphs and tables were used if available and approximate 
equations were written for those not availably. Graphs of the approxi- 
mations are shown in the appendix and test data should be checked for 
applicability of the approximations. 

Mass Rate of Air Flow . The mass rate of air flow, m, through an 
ASME square edged orifice is given in the Power Test Code {7] as: 



where 



m = 359 C F F^Y V^PqPo ’ 

C = coefficient of discharge 

F = velocity of approach factor 

d = orifice diameter, in. 

F^ = orifice thermal expansion factor 

Y = net expansion facto): 

APq = pressure differential across the orifice, in. H 2 O 

Po * specific weight of the air upstream of the orifice, 
lb„/ft3 
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The ratio of the orifice diameter to the duct diameter is: 



e = ^ 

D 

where D = duct diameter 

= 13.875 in. 

and the velocity of approach factor is: 



F = 



I 



The absolute duct pressure upstream of the orifice is given by: 



Pq = PbC0.490) - Po (0.0361) psic3 



where 



= barometric pressure, in. Hg 
Pq = duct pressure upstream of the orifice, in. H^O 



IS : 



(duct pressure is below atmospheric pressure) 

A dimensionless ratio of orifice pressure drop to orifice pressure 

APo (0.0361) 



X = 



Po 



and the net expansion factor is then: 



Y = 1 



,0 - ^.^yo.41 -P 0.3 5 



The approximate equation for the orifice thermal expansion factor 
for type 304 stainless steel gives: 

= 1.0+ (^O ^ 0 ^10 ) 

The coefficient of discharge is an unknown in the equation but may 
be quickly iterated using an equation suggested by Murdock [6]. The 



equation is dependant upon the Reynolds number in the orifice, ^nd is 



c - + AC 
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The coefficients 


Cq and AC 


for applicable 8*s 


from Murdock 


suggested first 


iteration C 


are: 




8 


Co 


AC 


C 


.15 


0.59446 


0.00945 


0.5975 


.25 


0.59483 


0.01037 


0.5966 


.45 


0.59863 


0.01543 


0.6014 


.75 


0.60480 


0.05448 


0.6128 



The density of the air upstream of the orifice may Ije found from 
the equation of state of a perfect gas modified by a humidity correction 
factor which gives: 



“ 



where 



14 4 fp Xpn 
53»3 4 Tq 

duct pressure upstream of the orifice, psia 



Tq = temperature of air upstream of the orifice, 

= humidity correction factor for density as given by 
Kays and London [A] 



0 



(I + 1.607 H') 



4n4 



H 



= humidity ratio, Ibs^j^ 

The dynamic viscosity, Mo, for air at one atmosphere pressure is 
given in reference (1) as: 



Mo - 



3/2 

- 0« 00 3 527 Tok 



Tok + 110.4 



where 



T^^ = temperature of air in duct, K 
The orifice Reynolds number may be found from: 



K, _ 15.2 8 m 
^ Red yWo d 
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Heat Transfer Calculations, The unit conductance for thermal con- 



vection heat transfer as given by el^uation one in the theory section is: 

1 



where 






j_ 

% 






1 



U 






4 - 



'“W 



12 



■>}o 

u 

^a 






"^wa 






= unit conductance for thermal convection heat 
transfer to the air, BTU/(hr ft^ 

= total surface effectiveness for the air side 

= overall unit conductance, BTU/ (hr ft^ ^F) 

= total heat transfer area on the air side, ft^ 

= total surface effectiveness for the steam side 

2 

= total heat transfer area on the steam side, ft 

= unit conduction for thermal convection on the 
steam side, BTU/(hr ft^ OF) 

= thickness of the wall separting the steam side 
from the air side, in. 

2 

= area of the wall on the air side, ft 
= thermal conductivity of the wall, BTU/ (hr ft^ 



The film coefficient for convection on the steam side, h^ , is unknown 
and must be estimated. Any reasonable estimate is considered adequate 
since a 50 percent error in the estimate results in about two percent 
error in h^. A value of hg equal to 2000 BTU/ (hr ft^ ^F) has been used. 

The specific heat of dry air may be considered constant and is assumed 

“ '^pda = 0.241 BTu/(lb„ 'f) 



The humidity correction factor for the specific heat as given by Kays 
and London [4] is: 

« I I > 91 5 H 
^c“ I + H 

where H = humidity ratio, lb H20/lbjjj dry air 
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I .r 




and the corrected specific heat is then: 

Cp = Xq Bxy/tibpr, f) 

The bulk average temperature in the test core is: 

■bi -t- to <?|- 

tb= ' 2" ^ 

where ti = air temperature upstream of the test core, op 
t2 = air temperature downstream of the test core, 
and the film average temperature is then: 

tf= 

^ 2 



where tg = the saturation temperature corresponding to the 
arithmetic average steam pressure in the core, 

The number of heat Transfer Units may be found from the temperature 

NTU relation as given in the theory section, yielding: 

and the overall .unit conductance is then: 



u 



- BTU/Chr "F) 

Aa 



For most plate-fin type geometries the equation for the temperature 
effectiveness of a straight fin with coi[istant conduction cross section 
may be used to a good approximation. [3] . The fin temperature effective- 
ness for the steam side, Hfg, is: 

tanh 



where 



Hg = effective fin length for the steam side and is 
equal to one-half the total fin height, in. 
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and hg = 2000 BTU/(hr ft^ °F) as previously estimated 

k^g = thermal conductance for the fin on the steam 
side, BTU/(hr ft^ °F/f t) . 

6g = fin thickness on the steam side, in. 

The overall temperature effectiveness for the steam side, rig, is then: 

^Wfi ~t" ^fs ^fs 

" As 

2 

where A^g = wall area on the steam side, ft 

2 

= fin area on the steam side, ft 
Ag = total heat transfer area on the steam side, ft^ 



The fin temperature effectiveness for the air side, "Hfa , is: 

tanh 

%a - 

where = effective fin length for the air side and is equal 

to one-half the total fin height, in. 

"1 



ma = 



24 ha 



ft 



^1 



and k£g = thermal conductance for the fin on the air side, 

BTU/(hr ft2 op/f t) . 

6^ = fin thickness on the air side, in. 

h^ = film coefficient for convection on the air side. 

It may be taken equal to U for the first iteration 
BTU/(hr ft2 op). 

The overall temperature effectiveness for the air side, Pq, is given by: 

^wa ^fa ^fa 
Aa 



'no = 



where 



wa 



= wall area on the air side, ft^ 



^fa = fin area on the air side, ft^ 



A^ = total heat transfer area on the air side, ft^ 
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The air mass velocity in the core, G, is given by: 



= A IbAfhr ft^) 

Ar 



where = the free flow area of the core on the air side, ft^ 

Stanton number in the core can now be evaluated as: 

Nst, = 7^ 

G Cp 

The dynamic viscosity evaluated at the film average temperature and 
one atmosphere pressure as given by Hilsenrath [1] is: 

// - 0.003527 'rk Ib^/fhr ft) 

Tfk+ 110.4 ’ 



where 



T^^ = film average temperature, 



The dynamic viscosity evaluated at the bulk average temperature is: 

3/2 

0.005527 Tbk ibpn/(hr ft) 

Tpk + 110.4 ' ' 



where T^^ = bulk average temperature, °K 

Ther thermal conductivity of air at the film average temperature 
and one atmosphere pressure as set forth by Hilsenrath [1] is: 




0.001529 “Tfk 



1/2 



24 5.4 



TfkXIO 






fk 



B TU 

hr °l/ft 



Where Tfk = f ilm average temperature, 

The Reynolds numbers evaluated at film average temperature and bulk 
average temperature are evaluated by: 



and 



NRef 



4hh6 



N 



Reb 
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where hydraulic radius r^ is defined by Kays and London [3] to be based 



on the free flow area and takes the form 

AcL 



>h = 



ft 



where L = core flow length, ft 
Prantdl number in the core is: 

_ 



Npr = 



kf 



where the fluid properties are evaluated at the film average temperature 
in the test core. 

The Colburn-j heat transfer modulus is then: 

• zA> 

J E Ns^Np/ 

which concludes the heat transfer calculations. 

Friction Factor Calculations . The equation for the core friction 
factor as developed in equation two in the theory section may be written 
as : 



f = 



■“h Pm 


APc (4.3255 X 10^) 


Kc -(i + <r^) 


L 


g2 


e, J 



K,? -h 1 (T" 



ez 



+ 4 



where r^^ 

Pm 

L 

APc 

G 



core hydraulic radius, ft 

mean density in the core Ib^/ft^ 

core flow length, ft 

pressure drop across core, in. H 2 O 

core mass velocity, Ib^/hr ft^ 
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cr 

Ke 

f2 

fd 



contraction coefficient 

the ratio of free flow area to frontal-area 
density of the air upstream of core, Ib^ft^ 
expansion coefficient 

density of the air downstream of core, Ibj^/ft^ 

friction factor for the duct downstream of the 
core 



1^ = length of duct from the core to the downstream 

pressure tap, in. 

^d ” hydraulic diameter of the duct downstream of 
core, in. 

The duct friction factor, fd, is assumed constant at .015, a value 



corresponding^ to a nearly smooth pipe for a Reyholds number range of 
100,000 to 500,000 as may be expected in the downstream duct. 



The mean density in the case of heating is given by: 




2 








where 



^2 

Ts 



Pi 



Pi 

Xm 

P2 



= upstream temperature, 

= downstream temperature, °R 

= Saturation temperature corresponding to the average 
steam pressure in the core ^R 

= upstream density 

= Q>0973I8 P| 

T| 

= upstream pressure, in. H 2 O absolute 
= humidity correction factor as previously given 
= downstream density 

^ 0.09318 Xm P2 Ib^^/ftS 
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downstream pressure, in. H 2 O absolute 



In the isothermal case, the mean density is equal to the arithmetic 



average density and since the temperature remains constant 

O _ e - 0-097318 X^(P| t P;?) Ib„/ft^ 



2T, 



m 



The equations for the coi^traction coefficient, K^, and the expansion 
coefficient, are developed by Kays [2] and are given by 



Kr = 



- 20c+ Oc 



(2Xd-0 



and 



c; 



Ke = I - 2Kd cr -H 0-‘ 



where = velocity-distribution coefficient 

<r = Ac/Afr 

= jet contraction ratio 
The jet contraction ratio is approximated as 

Cc = 0.611 -HO. 45 a- 4- 0.344 

The velocity-distribution coefficient Kj is dependent upon both core 
geometry and film Reynolds, • The usual geometries used in plate 

fin type construction are: 



square 

triangular 




For laminar flow, less than 2000, Kays [2] gives the values 

square Kj ” Kj =1.39 

U Ug 

triangular K^ = K^^ - 1.43 

' > L 

and for turbulent flow 
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square 



Kh = Kh = 1 + 1.17 (Kd, - 1) 



triangular Kd = Kd^. = 1 + 1.29 (Kd^ - 1) 

where the velocity-distribution coefficient for circular tubes, Kd^, is 
given by: 

Kcj^ ^ U09068 (4frr,) -f 0»0584 I ‘ 

- 0.2 



and 



— 0*049 



If the boundary layer is continually tripped, as in the case of louvered 
fin or off-set fin construction, film Reynolds number should be taken 
equal to infinity which yfelds Kd = 1.0 in all cases. 

Heat Balance Calculations . The heat transfer rate to the air, 

Q air, is: 



Q air = m Cp (t 2 “ ti) BTU/hr 

The total heat given up by the steam, Q steam, is equal to the heat 

given up by the excess steam as it passes through the core plus the heat 

given up by the condensate in the core. 

Q steam = s (hi h2g) + (hi - h 2 f) BTU/hr 

s = mass rate of flow of excess steam, Ibm/hr 

Wq = mass rate of flow of the condensate, Ib^j/hr 

= enthalpy of the steam at the top header pressure and 
temperature BTU/lb^ 

^2g” enthalpy of saturated steam at the bottom header 



where 



pressure, BTU/lb^ 



m 



h2f= enthalpy of saturated liquid at the bottom header 
pressure, BTU/lb 



m 



The mass rate of flow of the excess steam, s, as measured by a square 
edged, standard ASME orifice with flange pressure taps is found from the 
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applicable equation given in the ASME Power Test Code [7}: 



where 



I = 359 C 5 Fs d| F/i.,Ys ib„/ hr 

Cg = coefficient of discharge 

Fg = velocity of approach factor 

= orifice diameter, in.. 

= orifice thermal expansion factor 
= net expansion factor 



s 

^As 

Ys 

APs 

Vs 



= pressure differential across the steam orifice, 
in. H 2 O 

= specific volume of the steam passing through the 
orifice, ft^/lb^ 



The velocity of approach factor, Fg, is; 






where 



and 



Dg = steam pipe diameter 
= 1.25 in. 



The orifice thermal expansion factor is taken as constant corres- 
ponding to a steam temperature of 220 op. por type 304 stainless steel 



Fas = 1-003 



The net expansion factor for steam, Yg, is: 

Ys=|-(0.41 -h 0,35^#)('><s/l#3) 

where ' ' 

” APs / Ps 

The upstream steam pressure, Pg, is taken as the saturation pressure 
corresponding to the upstream temperature and is given by the approximate 
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equation 



Ps = 16.5 -f (tso -2l8Xo.376) psta 



The specific volume is given by the approximate equation: 

V 5 = Z4#0I - 0.4(tso *”218) fty lb 

The coefficient of discharge, C^, is iterated by the same method 
that was used for the air orifice. The values of and AC for applicable 
values of 3s are: 



.56 .700 in. 

.71 .890 in. 

and then 

C3 — Cos^~ 

The steam orifice Reynolds number, 

^ 0.004244 S 
•^Resd ds 

and the dynamic viscosity is given by the approximate equation: 

^5=8.25 X IO*^-l-(tso“200)(l.4 5 X 10®) 

The enthalpy of the superheated steam in the top header is approxi- 
mated by two equations, the first giving the variation with temperature 
at a constant pressure, Pgj^ = 18.0 psia and the second equation the 
variation with pressure. These equations are: 

hg^ = 1158.0 + (tg| - 230.0) (.5) BTU/lbjjj 

where h^j^ = enthalpy at 18.0 psia- 

tgl = top header steam temperatrue, 



ft sec 



~ *^os 



AC = AC, 



.6085 

.6090 



.025 

.041 



10 



N 



Resd 



r 



^Resd 3 



IS : 



35 



9 







thence 



hgi = - (Pgi - 18.0) (.35) BTU/lbm 

where Pg]^ = steam pressure in top header, psia 

The enthalpy of the saturated steam in the bottom header is found 
from the approximate equation: 

hg2 = 1155.3 + (Pg2 - 19.0) (.934) BTU/lb^ 

where Pg 2 = stream pressure in bottom header, psia 

and the enthalpy of the condensate in the bottom header is: 
h^ = 193.4 + (Pg2 - 19.0) (2.64) BTU/lb^ 

The mass rate of flow of the condensate, w^, is found by timing 
the condensate collected. 

5. Evaluation 

A steady state steam to air test was conducted on a compact heat 
exchanger manufactured by Hamilton Standard Corporation. The test core 
had the approximate dimensions of a six inch cube. The construction was 
the plane plate-fin type with flow passages between fins being approxi- 
mately square. One side of the core was “straight , ** meaning that the flow 
passage was straight for the length of the exchanger. The other side was 
“ruffled," meaning that the flow passage was approximately a cosine curve. 
The plate spacing was 0.156 inches and the fin spacing was 18 fins per 
inch for both sides. All plate and fin material was 0.004 inch stainless 
steel. A set of “standard" curves for Colburn-j and friction factor vs. 
Reynolds number were supplied with the test core. The “standard" curves 
were for plate and fin spacing similar to the core tested. 

The test was performed using the straight side of the core for air 
and the ruffled side for steam. The various comppnents of the facility, 
the air system, steam system and instrumentation were checked out and 
proved to be very satisfactory. There was a problem in defining flow and 
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heat transfer areas for the exchanger because structural stiffeners had 
been placed at the end of each plate-fin passage. The stiffeners were 
placed at an angle across the flow path and thereby restricted the 
flow in the peripheral flow passages. This restriction would tend to 
increase with increasing flow rates. The test results based upon the 
minimum free flow area, which neglects the peripheral flow passages 
affected by the stiffeners, was plotted on the standard curves as 
shown in Fig. 6. 

The next test should be on a heat exchanger core with previously 
determined characteristics and with given or definable geometry. 
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Figure 6, Colburn “3 and friction factor 
vs, test core Reynolds number. 
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APPENDIX I 



OPERATING PROCEDURE 

General . The preliminary check of instrumentation for both the 
hot core and cold core tests is the same. 

The multipoint recorder should be allowed to warm up two hours 
prior to calibration. To calibrate the instrument, set the print 
mechanism to ”Hold On** indicating any one of the points not used for a 
thermocouple. This open circuits the thermocouple measuring circuit. 
Place a D.C. millivolt source across the terminals at the back of the 
recorder and set the source to one millivolt. Adjust the recorder zero 
so that one millivolt is indicated. Set the source to nine millivolts 
and adjust the span so that nine millivolts are indicated. Repeat the 
zero and span adjustments until the recorder is calibrated. 

All manometers should be adjusted for zero with the isolation valves 
open prior to lighting off the air or steam systems. The two steam 
header mercury manometers must be zeroed at about three inches due to the 
head of water between the manometers and the water pots. The water lines 
between the mercury manometers and the water pots are susceptible to the 
formation of air pockets which will cause erroneous readings. 
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Cold Core Test> 



In order to assure an even distribution of data 



points, the approximate pressure differential across the air orifice 
should be calculated for each run. Calculate the mass rate of air 
flow at a core Reynolds number of 500 and 10,000 and plot these two 
points on a log log plot. Divide the straight line connecting the points 
into equal increments, the number depending upon the desired number of 
runs. The orifice pressure differential may then be read from Fig. 7 
for each mass rate of air flow. 

Open all manometer valves and adjust for zero. Close all manometer 
valves prior to lighting off. Set the multipoint recorder to print the 
thermocouples in use by raising the corresponding buttons on ”Selecto- 
Print.” Turn the chart motor on and select the desired speed. Close 
the compressor blast gate (discharge valve) and energize the compressor 
motor. Open the blast gate to desired setting. Adjust the sliding 
plate valve for the desired air orifice pressure differential and observe 
instrumentation until steady state obtained. 

Record the following data for each run: 

1. Atmospheric pressure (0.001 in, Hg) 

2. Atmospheric dry bulb temperature (0.1 ^F) 

3. Atmospheric wet bulb temperature (0.1 °F) 

4. Air orifice size 

5. Orifice upstream pressure (0.01 in. H^O) 

6. Orifice pressure differential (0.001 iru H 2 O for APo <3”, and 
0.01 in. H 2 O for AP^ >3”) 

7. Core upstream pressure (0.01 in. H 2 O) 

8. Core pressure differential (0.001 in, H 2 O for APc <3**, and 0.01 
in, H 2 O for APc >3**). 
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9* Orifice air temperature (0.01 millivolt) 

10. Core upstream air temperature (0.01 millivolt) 

Items five through ten are recorded three times in the same sequence 
and averaged on the data sheet. 

Hot Core Test . The instrumentation checks and the procedures for 
lighting off the air system are the same as for the Cold Core Test. 

The steam system is activated by the following procedures. Obtain 
100 psig low pressure air to the reducer central valves. Turn on the 
tap water supply for the condensate cooler and the desuperheater. Crack 
the desuperheater supply valve. Open the by-pass valves for the steam 
reducers and crack the steam stop valve. Allow the system to warm up 
slowly. When warm-up is completed close the by-pass valves and open the 
steam stop. Charge the first steam reducer to about 30 psig and the 
second to about six psig. Obtain a steam state in the top steam header 
of about six psig with five to ten degrees superheat by adjusting the 
first steam reducer and the desuperheater. The excess steam flow rate 
should be at least five times the condensate flow rate. Two steam 
orifice plates are provided; the larger to be used for high air flow rates. 
Condensate should be collected for at least five minutes at high air flow 
rates and for at least ten minutes at low air flow rates. 

Steady state must be achieved before each run. This will require 
about one hour prior to the first run and 15 minutes between runs. Steady 
downstream core temperatures are the best indication of steady state. 

The following data is recorded for each run: 

1. Atmospheric pressure (0.001 in. Hg) 

2. Atmospheric dry bulb temperature (0.1 °F) 
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3. Atmospheric wet bulb temperature (0.1 op) 

4. Air orifice size 

5. Steam orifice size 

6. Top steam header pressure (0.01 in. Hg) 

7. Bottom steam header pressure (0.01 in. Hg) 

8. Pressure differential across the steam orifice (0.01 in. Hg 
under H 2 O) 

9. Core upstream pressure (0.01 in. H 2 O) 

10. Core pressure differential (0.001 in. H 2 O for aPc< 3”, and 0.01 
in, H 2 O for APc> 3**) 

11. Orifice pressure differential (0.001 in. H 2 O for APo< 3”, and 
0.01 in.H20 for APo >3”) 

12. Orifice upstream pressure (0.01 in. H 2 O) 

13. Core upstream air temperature (0.01 millivolt) 

14. Core downstream air temperature (0.01 millivolt) 

15. Orifice air temperature (0.01 millivolt) 

16. Top steam header temperature (0.01 millivolt) 

17. Bottom steam header temperature (0.01 millivolt) 

18. Steam orifice temperature (0.01 millivolt) 

19. Condensate weight (0.1 Ih) 

20. Condensate collection time (1. sec^ 

Record items 6 through 18 three times in the same sequence and 
average on the data sheet. 

The humidity ratio, H, may be found from a standard psychometric 
chart and convert temperatures to (0.1). 
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Figure 7* Mass rate of air flow vs« 

orifice pressure differential. 
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APPENDIX II 



DIGITAL COMPUTER PROGRAM FOR DATA REDUCTION 
The program is designed to reduce data from both the hof core test 
and the cold core test. The program written in FORTRAN 60 and is shown 
on pages 32 through 57 and the output for the hot core test is shown on 
page 58 and that of the cold core test on page 59 . The variables used 
in the program are defined in the program glossary. Variables not appear- 
ing in the glossary are internal to the program and are used to change 
dimension or to define special groupings to simplify programing. 

The index used to reduce data from the hot core test or the cold 
core test is given by: 



Program Index 
IHORC 

Hot core test 1 

Cold core test 2 

The program is indexed to accept the various fin geometries normally 
used in plate fin construction and the indices; are as follows: 

. Fin Geometry Index 
ISQTRIN 

Square fin 1 

Triangular fin 2 

Louvered or off-set fin 3 

(Nj^^ = infinity) 

"The indices used for the fluid metering orifices are: 

Air Orifice Air orifice Index 

BETA d (in.) IB 
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.75 



10.406 



4 



.45 



.25 



.15 



Steam Orifice 
BS 



.71 



.56 



6.244 3 

3.468 2 

2.081 1 

Steam Orifice Index 
ds (in.) IBS 

0.890 2 

0.700 1 



The indices are all entered as fixed point constants. The number 
of runs, NOR, and the run number, NR, are entered as fixed point con- 
stants. All fixed point constants must be right justified in their 
data fields. 

The constants for a given core such as areas, lengths, etc. are 
entered as input data only once and are the first three data cards 
beginning with AC and ending with DD. The fixed point constants for a 
particular series of runs, NOR, IHORC, ISQTRIN are entered on the fourth 
data card. The fifth through the seventh data cards are for run one, 
eighth through tenth for run two and so on. The first data card for a 
run contains the run fixed point constants, NR, IB, IBS. The second and 
third data cards for a run contain the run data beginning with H and 
ending with TC. When entering data for the cold core test, only the 
applicable data must be entered. All other data fields may be left blank. 

The dimensions of input data must be those given in the program 
glossary and are the same as the raw data. 
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PROGRAM GLOSSARY 



Program 

S3mibol 


Program 

Unit 


Nomenclature 

S}mibol 


AA 


ft^ 


^a 


AC 


ft2 




ALFA 






AFA 


ft2 


^fa 


AFR 


ft2 


Afr 


AFS 


ft2 


Afs 


AKFM 


BTU/(hr ft2 °F/ft) 


Kfm 


AS 


ft2 


As 


AWA 


ft2 


Awa 


AWS 


ft2 


Aws 


BETA 




3 


BS 






C 




C 


CC 




Cc 


CO 




Cq 


COLJ 




j 


COS 




^os 


CP 


BTU/(lbjn °F) 


c 

p 


CRKD 




dc 


CS 




Cs 


DC 




AC 


DCS 






DD 


in. 


Dd 
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DL 


in. 




DPC 


in. H 2 O 


> 

0 


DPO 


in. H 2 O 


AP 

0 


DPS 


in. Hg under water 


> 

CO 


DSO 


in. 




ERROR 




Error 


FA 




Fa 


FAS 




Fas 


FB 




F 


FBS 




Fs 


FC 




f 


FD 




fd 


FLA 


in. 




FLL 


ft 


L 


FLS 


in. 


CO 


FKA 


BTU/(hr ft^oF/ft) 


kfa 


FKS 


BTU/(hr ft2op/ft) 


CO 


FM 




Fm 


G 


lbs/ (hr ft^) 


G 


H 


lb /lb 


H 


HA 


BTU/(hr ft2op) 




HG 


BTU/lb„ 

m 


h 

c 


HS 


BTU/(hr ft^op) 


hs 


HSl 


BTU/lb 

m 


hsl 


HS18 


BTU/lb^j 


^sl 


HS2 


BTU/lbjjj 


hs2 
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IB 


Index for air orifice diameter 


IBS 


Index for steam orifice diameter 


IHORC 


Index for Hot or Cold core test 


ISQTRIN 


Index for fin geometry, square, triangular 
or infinite Reynolds number 


NOR 


Number of runs 


NR 


Run number 


OD 


in. d 


PI 


in. H 2 O Pj^ 


P2 


in. H 2 O P 2 


PB 


in. Hg 


PO 

PR23 


in. H 2 O Pq 

N 2/3 
Pr 


PSI 


in. Hg Pgi 


PS2 


in. Hg Pg2 


PSO 


Psia Pgo 


Q AIR 


BTU/hr Q air 


QSTM 


BTU/hr Qsteam 


REB 


%eb 


REFM 


%ef 


REO 


\ed 


RESO 


^Resd 


RH 


ft r^ 


ROl 

R02 


lb' P]^ 

lb, P2 


ROAVE 


m/^^ Pave 
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ROM 



ROO 

SDOT 

ST 

SQKD 

T1 

T2 

TB 

TC 

TFA 

TFM 

TFS 

TO 

TS 

TSl 

TS2 

TSO 

TRKD 

TW 

U 

UB 

UC 

UFM 

UO 

USO 



lb „/f 


Pm 


lb 


Po 


lb ^/hr 


s 




Nst 




Kd 




s 


Op 




Op 


^2 


Op 


tb 


min 


time of condensate 
collection 


in. 


<Sa 


Op 


tf 


in 


<5s 


Op 


to 


Op 


ts 


Op 


^sl 


Op 


ts2 


op 


^SO 






in. 


1 

w 



BTU/(hr f2°F) U 

lb' In/ (hr ft) 

lb ni/(hr ft) yU in core for cold core 

test 

lb (hr ft) 
lb^m/(hr ft) 
lb m/(sec ft) 
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vs 



wc 

WCH 

WDOT 

WK 

X 

XC 

SKD 

XM 

XMA 

XMS 

XNTU 

XS 

Y 

YKE 

YNFA 

YNFS 

YNO 

YNS 

YS 

Z'KC 



ft3/lbm 


Ys 


Ib.m 


Weight of ' condensate 
collected 


lb ^/hr 


Wc 


lb ta/hr 


# 

m 


BTU/(hr ft^OF/ft) 






X 




X 




c 




^d 




Xm 


ffl 


nia 


ft-1 






NTU 




Xs 



Y 

s 

Ys 
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PROGRAM bSHEAT 

OREAD 1,AC,AFR,AA,AWA,AFA,AS/AWS*AFS,FLL»TFS/TFA»TW, 

ifla,fls,fka»fks,wk,fd»dl,dd 

1 FORMAT(8F10.0) 

READ 2, nor, I hORC, ISOTRIN 

2 FORMAT(3I10) 

GO TO (3,5) , IHORC 

3 PRINT 4 

40FOKMAT(1H1 25X,32H HEAT TRANSFER AND FRICTION DATA,/// 

1106H RUN MOOT G HUMID T1 PI T2 DPC T^ 

2H BTU/ COLJ F FACT NRF NR B ERROR,// 

3107H IB/hr lb/hr FT2 L8/LB F PSIA F PSI F H 

4R FT2 F PERCENT,//) 

GO TO 7 

5 PRINT 6 

60FORMAT(1H1 25X,25H ISOTHERMAL FRICTION DATA,/// 

160H RUN MDOT G HUMID Tl PI DPC F FACT NR,// 

246H LB/HR LB/HR FT2 LB/LB F PSIA PSI,//) 

7 READ 8, NR, IB, IBS 

8 FORMAT (3110) 

OREAD 9, H, OD, DSO , BET A , 8S , TO, T 1 , t 2 , TSl , TS2 , TSO , 

1PB,PS1,PS2, DPS,P1,DPC,PO,DPO, WC,TC 

9 FORMAT (8F10.0) 

PS1A= . 490*(PB-^PS1) 

PS2a=.490*(PB-^PS2) 

PSAVE=(PS1A*PS2A)/2. 

TS=220.f(PSAVE-l7.1)*2.66 
p2 = Pl-*-DPC 
T0R=T0>459.7 
T0K=10R*5./9. 

UO= . 003527* TOK^^l ,5/(TOK+110 .4) 

GO TO (20,21,22,23 ),IB 

20 CO=. 59446 
DC= .00945 
Cl=.5975 
OD=10 • 406 
BETA= ,75 
GO TO 25 

21 CO=. 59483 
DC=. 01037 
Cl= . 5966 
0D=6.244 
BETA=. 45 
GO TO 25 

22 CO=. 59868 
DC= .01543 
C1=.6014 
0D=3 . 468 
8ETA=.25 
GO TO 25 

23 CO=. 60480 
DC=. 05448 
Cl=.6128 
0D=2 . 081 
BETA=,15 

25 C=Cl 

B4=8ETA**4 

POA = PB*, 490’'PO*,0361 
X=(DPO*,0361)/POA 

Y=l. -( ( ( .41+.35 *b4)*X)/1. 4) > . . 

XM=(1.*H)/(1.*1,607*H) 
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ROO=POA*144.*XM/(5i.3*TOR) 
IF(T0-68. )12,12,11 

11 FA=l.+(T0-68. )*1.85E-5 
GO TO 13 

12 FA=1. 

13 F8=1,/(1.-B4)**.5 



DO 26 I=l#3 

WDOT=359.*C*FB*FA*Y*OD*OD*SQRTF(DPO*ROO) 

RE0=15,28 *WDOT/(UO*OD) 

C = CO-*-DC*SQRTF( 10 0 0 0 ./R60) 

26 CONTINUE 

GO TO (29,40 )# IHORC 
29 TB=(Tl+T2)/2. 

TFM=(TB+TS)/2. 

TBK=(T6+459.7)*5./9. 

TFMK=(TFMf459.7)*5./9. 

UB=.003527*TBK**1,5/(TBK+110.4) 

UFM= . 0 03527*TFHK**1 .5/ ( TFMK-^110 . 4 ) 

AKFM= .0 01529*TFMK**. 5/(1. +245.4/ (TFMK*10 .**(12. /TFMK) ) ) 

RH= AC*FLL/AA 
G=WD0T/AC 

refm=''.*hh*g/ufm 

rEB=4 . *RH*G7UB 
XC=(1.+1.915*H)/(1.+H) 
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CP= XC*,241 

XNfU =LOGF( (TS-T1)/(TS-T2) ) 



u=cp*wdot*xntu/aa 
HS=2000 , 

XMS=SQRTF(24. *hS/(FKS*TFS) ) 

YNFS=(TANH(XMS*FLS/12.))/(XMS*FLS/12.) 
YNS=( AWS+YNFS*AFS)/AS 
HA = U 
DO 30 1=1,3 



XMA=SQRTF(24.*HA/(FKA*TFA) ) 
YNFA=(TANH(XMA*FLA/12.))/(XMA*FLA/12.) 

YN0=(AWA + YNFA*AFA)/AA ^ HO VXXN 

HA = l./(YN0*((l./t))'-AA*(l./(YNS*AS*HS) + TW/(AWA*WK*12,)))) 



PR23= (UFM*CP/AKFM)**.667 



ST=HA/ (G*CP ) 
COLJ= ST*PR23 



40 ALFA= AC/AFR 

P1A= PB*13.571- PI 
P2a= PB*13.571- P2 
T1A= Tl+459.7 
T2A= T2+459.7 

R01=. 097318* XM*P1A/T1A 
R02= . 097318*XM*P2A/T2A 

GO TO (41,42) , IHORC 

41 R0M=.5*( (TlA*R01 + T2A*R02)/(TS + 459.7-( (T2A’T1A)/XNTU) ) ) 

GO TO 43 ^ . 

42 ROAVE = .097318*XM*(P1A+P2A)/(2,*T1A) 

R0M=R0AVE 

G=WD0T/AC 

rh=ac*fll/aa 

T1K=(Ti+459.7)*5./9. 

UC= .003527*TlK**1.5/(TlK+110,4) 

REC=4.*RH*G/UC 

REFM=REC 

43 GO TO (50,45,44) , ISQTRIN 

44 XKD=1. 

GO TO 60 

45 IF(REFM-2000. )46,47,47 
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46 TRKD=1.43 
GO TO 48 

47 FM =.049*REFM**(-,2) i : 

CRKD = 1 . 09068*4 , *FM+ . 05884*SQRTF( 4 . *FM) +1 . 

TRKD=1.+1.29*(CRKD-1. ) 

48 XKD=TRKD 
GO TO 60 

50 IFCREFM -2000 . )51,53,53 i 

51 SQKD = 1.39 

52 GO TO 56 ; 

53 FM =.049*REFM**(-.2) 

54 CRKD=1. 09068*4. *FM*. 05884*SQRTF(4,*FM)*1. ‘ 

55 SQKD = 1.*1.17*(CRKD-1. ) i: 

56 XKD=SQKD \ 

60 CC = . 611* . 045*ALFA * . 344* ( ALFA**5 . 7 ) i 

YKE = 1.-2.*XKD*ALFA*ALFA**2. j. 

ZKC=(l.-2.*CC*(CC*CC)*((2.*XKD)-l.) )/(CC*CC) | 

A=(ZKC- (1.*ALFA*ALFA) )/R01 • 

B=(YKE+1.+ALFA*ALFA*(1.*(4.*FD*DL/D0)))/R02 

FC= (RH*R0M/FLL)*( (4.3255E*9*DPC/(G*G) )-A-B) ■ 

PIAA = P1A*.0361 ; 

DPCA=DPC* . 0361 

GO TO (70,104) , IHORC I 

70 GO TO (75,76), IBS I 

75 COS=.6090 
DCS=.0410 
DSO= .89 
BS=.71 

GO TO 77 

76 COS=.6085 
DCS=.0250 
DSO=.70 

BS=.56 ; 

77 CS=.61 

BS4=BS**4, f 

FBS=1./(1.-BS4)**.5 

FAS=1.003 

DPSO= . 454*DPS c 

PS0 = 17 . l+( TSO-220 . ),*.376 
XS=DPSO/PSO 

YS=1.'( .41*.35*BS4)*(XS/1.3) 

DPSW=DPSO/. 0361 
VS=24.01-(TSO-218. )*.4 
US0=8 . 25E-6+( TSO-200 . )*1 . 45E-8 
DO 80 IS=1,3 

SD0T=359.*CS*FSB*FAS*YS*DS0*DS0*SQRTF(DPSW/VS) 

RESO=.004244*SDOT/(USO*DSO) 

CS=COS*DCS*SQRTF ( 10000 ./RESO) 

80 CONTINUE 

HS18=1158.0*(TS1-230. )*.5 
HS1=HS18-(PS1-18. 0 )* .35 
HS2=1155.3*(PS2A-19. )*( .934) 

HC=193.4*(PS2A-19. )*2.64 
WCH=WC*60 ./TC 
QAIR=WD0T*CP*(T2-T1) 

QSTM=WCH* (HSl-HC )*SD0T*(HS1"HS2) 

ERROR=(QAIR-QSTM)*100 ./QAIR 

99 PRINT iOO,NR,WDOT,G,H,T1,P1AA,T2,DPCA,TS,HA,COL.J,FC,REFM,REB,ERROR 
1000FORMAT(I3,2X,F6.1,2X,F7.1,2X,F5.4,2X,F4.1,2X,F5,2,2X,F5.1,2X,F4,2, 
12x,F5.1,2x,F4.1,3x,F6.5,2x,F6.5,2x,F6.1,2X,F6,1,2X,F5.1,//) 

GO TO 110 - . 

'.04 PRINT 105,NR,WDOT,Q,H,T1,P1AA,DPCA,FC,REFM 

' SA 



1050FORMAT( 13,2 X,F6, 1,2 X,F7.1 

12X,F6*1,/ 

110 IF(NR-N0R)7,112,112 
112 CONTINUE 
END 



2X,F5.4,2X,F4.1,2X,F5.2,2X,F4,2,2X,F6,5» 
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APPENDIX III 



GRAPHS FOR APPROXIMATE EQUATIONS 



Graphs for the approximate equations 
used in the experimental method and 
equations section are given on the 
following pages. 
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Figure 8. Thermal expansion factor 
vs. temperature. 
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Figure 10. Saturated steam pressure 
vs . temp e ratur e • 



63 



/b/ft sec 




Figure 11. .Saturated steam dynamic 

viscosity vs. temperature. 
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Figure 12. Saturated steam specific 
volume V8. temperature. 
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Figure 15 • Saturated steam enthalpy 
vs. pressure. 
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Figure lA. Superheated steam enthalpy 
vs. temperature. 
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Figure I5. Saturated liquid enthalpy 
vs. pressure. 
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